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Abstract 

Rotational motion of extended celestial bodies is discussed in the framework 
of the Parametrized Post-Newtonian (PPN) formalism with the two param- 
eters 7 and (3. A local PPN reference system of a massive extended body 
being a member of a system of N massive extended bodies is constructed. 
In the local PPN reference system the external gravitational field manifests 
itself only in the form of tidal potentials. Rotational equations of motion, 
which are then derived in the local reference system, reveal a special term in 
the torque analogous to the Nordtvedt effect in the translational equations of 
motion: it is proportional to 4/3 — 7 — 3, to the acceleration of the body rela- 
tive to the global PPN reference system and to some quantity characterizing 
the distribution of inertial gravitational energy within the body. This term is 
a direct consequence of the violation of the Strong Equivalence Principle in 
alternative theories of gravitation. 

PACS numbers: 04.25.-g, 04.25.Nx, 04.80.Cc, 95.10.Ce 
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I. INTRODUCTION 



In present paper we deal with rotational motion of a celestial body being a member of 
a system of N extended bodies in the PPN formalism. It is well known that rotational 
motion of non- isolated extended celestial bodies is a subtle point. In General Relativity 
this problem has been treated in the post-Newtonian approximation by Fock |l[ who used 
one global reference system to describe both translational motion of N extended rotating 
bodies and their rotational motion. The resulting equations of motion derived by Fock 
showed complicated couplings between translational motion of the bodies and their rotation. 
However, most of these couplings were only spurious coordinate effects which were later 
shown to vanish when using physically adequate local reference systems to describe the 
local dynamics (e.g., rotation) of each body of the system (see M for further comments and 
references). The idea of the local reference systems is to introduce a reference system for a 
body in which the influence of external masses is effaced as much as possible. In General 
Relativity the local reference systems have been discussed in great details by Brumberg and 
Kopeikin (see also ||) and Damour, Soffel and Xu ||[7|j2|. In particular it has been 
shown that the local reference systems in General Relativity have two properties: 

(A) the gravitational field of external bodies is represented in the form of tidal potentials 
being C(X 2 ), where X 1 are local spatial coordinates; 

(B) the internal gravitational field of the body coincides formally with the gravitational 
field of a corresponding isolated source provided that the tidal influence of the external 
matter is neglected. 

The first attempt to use a version of such a local reference system to study the rotational 
motion of an extended body has been undertaken by Voinov ||. Major progress to solve 
the problem in the Einsteinian post-Newtonian theory was achieved by Damour, Soffel and 
Xu || who used their DSX formalism aimed at constructing the local reference systems 
and derived the rotational equations of motion of each body of an iV-body system with full 
multipole structure. 

The aim of this paper is to derive rotational equations of motion of an arbitrarily-shaped 
extended body being a member of a system of N arbitrarily-shaped extended bodies in 
the framework of the PPN formalism with the two Eddington parameters f3 and 7. The 
PPN formalism is a phenomenological scheme giving in generic parametrized form the post- 
Newtonian metric tensor for an isolated material source in a wide class of metric theories of 
gravitation. Many aspects of testing General Relativity in the weak-field slow-motion regime 
are based upon the PPN formalism. On the other hand, the theory of local reference systems 
is also proved to be very important for physically meaningful modeling of observational 
data. Therefore, it is quite important to elaborate a theory of local reference systems in the 
framework of the PPN formalism (see also ||). 

In full detail our theory of local PPN reference systems will be discussed elsewhere. Here, 
we confine ourselves to a brief description of the theory and to one of its results: rotational 
equations of motion of a body relative to its local reference system. These equations reveal 
an interesting physical effect. They contain an additional torque analogous to the Nordtvedt 
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effect in the translational equations of motion: the torque is proportional to the Nordtvedt 
parameter 77 = 4/3 — 7 — 3, to the acceleration of the body relative to the global PPN 
reference system and finally to some integral over the volume of the body characterizing the 
distribution of the internal gravitational energy inside the body. 



II. GLOBAL PPN REFERENCE SYSTEM 

We begin with the global PPN metric tensor of a isolated N body system in the form 
(see, e.g., Q) 

2 2 
£00 = -1 + -J w(t, x) - - (5 w 2 (t, x) + C(c~ 5 ), 



2(1 + 7) 

goi = - 1 T R w%x)+0(c 



-5\ 



C 



where 



w = G J ^^O^^^' + ^Gj^ J<T(t,x.')\xL-x.'\d 3 x' + 0(c- A ), (2) 
w { = G f a l (t,x')- — l — dV + 0(eT 2 ), (3) 



x — X 



i (V 00 + 7 T fcfc + ir 00 (3 7 - 2/3 - 1) ™) + o( c - 4 ), a 1 = 



and T a/3 is the energy- momentum tensor in the global PPN reference system (t, x l ). We 
retain here only two PPN parameters: 7 and (3. This simplest version of the PPN formalism 
covers, nevertheless, most viable theories of gravitation. For some technical reasons we prefer 
to work not in the standard PPN gauge, but in a different one, the difference affecting only 
the coordinate time t 

*PPN = t-^X + 0(cT 5 ), (5) 

where tppN is coordinate time of the global PPN reference system in the standard PPN 
gauge and 
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-G J a(t,x') Ix-x'IdV. (6) 



For 7 = (5 = 1 the reference system (t, x l ) is harmonic. 

We have written the metric tensor of the global PPN reference system without specifying 
an explicit form of the energy-momentum tensor. This is no more than a formal way to 
specify the global PPN metric tensor. The metric coincides with the version of the 

PPN formalism described in [n| as well as with the version discussed in |TTJ , provided that 



only two parameters 7 and f3 are retained in both versions. Effectively we consider only 



3 



those theories of gravitation which produce the metric ([!])— (]4]) in the first post-Newtonian 
approximation, the energy-momentum tensor being left unspecified. 

We consider the material system to consist of N bodies which represent simply spatially 
bounded blobs of matter. Potentials w and w l are defined by (ID)-© as volume integrals 
over the whole space. Using this fact, we split the area of integration into the volume V of 
the body, for which we want to construct a local reference system, and the remaining part 
of space. Thus, we split w and w % into internal potentials (potentials of the body under 
consideration) and external ones (potentials due to the other bodies): 

w(t, x) = w E (t, x) + w(t, x), 

w%x) = «4(t,x) + x). (7) 



III. LOCAL PPN REFERENCE SYSTEM 

From a physical point of view any reference system covering a space-time region under 
consideration can be used to describe physical phenomena within that region. Some reference 
systems, however, offer a simpler mathematical description for physical laws and make the 
character of physical laws more obvious. It is well known that, for example, for a massless 
observer one can construct a local reference system where gravitational fields appear only 
in tidal form. Due to the works of Brumberg and Kopeikin |^ and Damour, Soffel and 
Xu we know that in the first post-Newtonian approximation of General Relativity it 

is possible to construct an analogue of such a local reference system for a massive extended 
body and that those local reference systems have two properties A and B listed in the 
Introduction. It is clear that the possibility to satisfy in General Relativity both properties 
simultaneously is closely related with the Equivalence Principle. It can be expected apriori 
and will be seen below that in alternative theories of gravitation where the Equivalence 
Principle is violated properties A and B cannot be satisfied simultaneously. 

We assume the metric tensor of the local PPN reference system (T, X a ) of a selected 
body to be of the form 

Goo = -1 + \ W(T, X) - - (5 W 2 (T, X) + 0(c" 5 ), 

G a = - 2(1 t 7) W a (T, X) + 0(c" 5 ), 
cr 

Gab = Sab (l + 1 7 W(T, X)) + 0(c" 4 ), (8) 
and the local gravitational potentials W and W a to admit a representation in the form 

W(T, X) = W E (T, X) + Q a (T) X a + W T (T, X) + 1 *(T, X), (9) 

W a (T, X) = W£(T, X) + i e abc C b (T) X c + W?(T, X). (10) 

Here, the local internal gravitational potentials W E and W E have the same functional form as 
their global counterparts u>e and w E , but all quantities should be taken in the local reference 
system, i.e. 
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W E = G J v S(r, X') j rf 3 X' + ^G— ^S(T,X')|X-X'|d 3 X' + 0(c- 4 ), (11) 
W| = G J v Z a (T,X') ^ d 3 X / + 0(c- 2 ), (12) 

e = 1 (V 00 + 7 r aa + (3 7 -2/3-1) + e>( c - 4 ), s a = -r 0a + e>( c - 2 ), (13) 

c 2 V cr / c 

where T a/3 is the energy-momentum tensor in the local reference system. The external po- 
tentials Wt and W% represent tidal field of the other bodies of the system and are assumed to 
be ~ C(X 2 ). It is also assumed that (Wt, W%) are functions of (w,Hf) and their derivatives 
as well as of the trajectory of the origin of the local reference system relative to the global 
one. Two arbitrary functions Q a {T) and C a (T) have a clear physical meaning which will be 
discussed below. Finally, the function \?(T, X) is some unknown function containing internal 
potentials of the central body which appear in the local PPN reference system in addition 
to W E and W E . Clearly the appearance of $ is related with a violation of the Equivalence 
Principle which make it impossible to satisfy simultaneously properties A and B formulated 
in the Introduction. By assuming that W^ and W% are ~ 0(X. 2 ) we assume that property 
A is satisfied. Therefore, property B is violated which results in the appearance of \E r . 

Then, the results of the Brumberg-Kopeikin and DSX formalisms (see, e.g., Theorems 
1 and 2 of ||) allow us to write the transformations between the global and local reference 
systems in the form 

T = t-\(A + + -AB + BV£j + B^r E r E + C(t, x)) + 0(cT 5 ), (14) 



9 V 1 "KtW 1 A 



X a = U + - 2 (o< 44 + & k r k E + 7 (r E a E r E - - a E r 2 E ) ) ) + 0{c~% (15) 



1 

^2 \ 2 ^ ^ ^ ' £j ' ' I ' Hj Hi ' Hi ^ 

where r E = x l — x E (t), x E (t) is the coordinates of the origin of the local reference system 
relative to the global one, and v E = dx E /dt and a l E = d 2 x E /dt 2 are its velocity and acceler- 
ation, respectively. The functions A(t), B(t), B { (t), B^(t) = B ji (t), D^{t) = D^(t), R^t) 
(being orthogonal matrix) and C(t, x) ~ 0(r E ) are some unknown functions. 
The transformation rule 

r)X^ r)X u 

^(t,x) = ^ r ^- r G„(T,X) (16) 

enables us then to derive or constrain the unknown functions from both the local PPN 
metric tensor (|8|)-(|i"3D and the transformations (n)-(15) 

A(t) = ~vl + w(x E ), (17) 



B(t) = ~4 + 2( 7 + 1) 4^(x E ) v 2 E w(x E ) + (P~1) ^ 2 (xe), (18) 

1 



B\t) = ~v 2 E v E + 2(1 + 7) W l (x E ) - (2 7 + 1) 4«J(x E ), (19) 
D ij = 7^'m7(x e ), (20) 
BV(t) = -4^3)Qa + (1 + 7)^ (ij) (x E ) - 7 4 i ^ i) (x E ) + J 7 <^#(x E ), (21) 



5 



c 2 R\ R a . = -(l + 7 ) e ljk R\C a 

- 2(1 + 7) W^(x E ) + (1 + 2 7 ) 4 Wj^Xe) + 4#j]<2a + 0(c~ 2 ), (22) 

C« = ( 7 -2)a*r|, (23) 

4(f) = x E (t*)) - R a 3 Q a {T) (> - 1 (vl S ij + (2 + 7 ) ™(x E ) ^ + ^44)) 

+ -^2(1 + 7) #(x E ) + ( 7 4 - 2(7 + P) w(x E )) W,i(x E ) - (2 7 + 1) 4^(xe) 

- 2(1 + 7) 4 TOj(x B ) - 44^(x E )) + 0(cT 4 ), (24) 

T = t* - -A{t*)+0{ C - A ). (25) 

For any function A(x E ) means A(t, x E (i)). Parentheses and brackets around a group indices 
signify, respectively, symmetric and antisymmetric parts of the corresponding expressions: 
AM = \ (A ij + A ji ) and A [ij] = \ (A ij - A ji ), etc. The moment of time t* which appears in 
(p4|) is defined by ( p5|) . The function C(t, x) is not fixed completely, but should only satisfy 
(p3|). It is clear, however, that the PPN equations of motion do not depend on C(t, x). 
Equations fl2~4])-(|2"5"D represent equations of motion of the origin of the local reference system 
relative to the global one. For Q a = the equations coincide with the geodesic equations in 
the external potentials w and w 1 . 

The external potentials and W% represent the tidal fields of external masses and are 
given by 

W T (T, X) = w(t, x) - u>(x E ) - W i (x E )4 

+i (-2(1 + 7 )4 (^(*,x) -TZJ*(x E ) -^(x E )4) + (1 +7)4^t 

+ (1 + 7 )# J '(x E ) r E 4 + \ 7^(xe) 4 + Q - - 7) (44f 

+ (1 - 2/5 - 2 7 ) Q a X a a E r E - 7 44 (x E ) 4 + ^7 4 R%^Q a 

+ Ac(r, X) + 2(1 - 0) (w(xb) + 44) W T ) + 0(c~ 4 ), (26) 

W?(T, X) = i2° 4 |^(t, x) - ^(x E ) - W^(x E )4 - <W T (T, X) 

+ ^j-^y (7 (444 - 5 44) - C,i(T, X)) } + O(c^). (27) 

Finally, the function ^ reads 

*(T, X) = -r) (w K {t, x) (W(x E ) + 44) - x E 4 (t, x) a E ) + 0(cT 2 ), (28) 

77 = Af3 — 7 — 3 being the Nordtvedt parameter, and x £ (^ x ) is defined by the integral 
(§) taken over the volume of the central body. The presence of this function reflects 
the violation of property B in the local reference system (T, X a ). Property A is satisfied, 
since both Wt and W% are ~ 0(K 2 ). One can show that by changing the transformations 
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(fU])-(l25D one can construct another version of the local reference system where property 
B is valid, but property A is violated. One can also show that it is impossible to satisfy 
both properties simultaneously. This fact is a direct consequence of the violation of the 
Equivalence Principle for rj ^ 0. 

The function Q a (T) characterizes the world line of the origin of the local reference sys- 
tem. Q a (T) represents the acceleration of the instantaneous locally inertial reference system 
(whose origin coincides with that of the local reference system at a given moment of time) 
expressed in the local reference system. E.g., for Q a = the origin moves along a geodesic 
in the external gravitational field. The value Q a can be also chosen is such a way that the 
multipole expansion of the internal gravitational field W-^ does not have a dipole compo- 
nent. The function C a (T) describes the spatial orientation of the local reference system with 
respect to the global one. Its relation to the orthogonal matrix is defined by (|22|). One 
can choose C a (T) so that R a i = 5" and the resulting local reference system does not rotate 
relative to the global one. Another possible choice is C a = resulting in a dynamically 
nonrotating local reference system and the orthogonal matrix R a i in the transformations 
of the spatial coordinates represents the well-known de Sitter, Lense-Thirring and Thomas 
precessions. 



IV. ROTATIONAL EQUATIONS OF MOTIONS 



The local reference system of an extended massive body described above allows us to 
derive the rotational equations of motion of the body in the framework of the PPN formalism. 
From the equation (here, Q is the determinant of the local metric tensor) 

e abc [ (-g)X b T c %d 3 X = 0, (29) 

which is valid due to the local equations of motion 

T a % = 0, (30) 

in analogy to General Relativity (see one can derive the rotational equations of motion 

for the body 

^ S a = inertial + + ^Nor + (31) 

where the PPN spin S a is defined by 
S a = e abc J v X b p c (T,X) d 3 X + 0(c" 4 ), (32) 



c 



57 - 1 „ A 1 ^ f VA (47 + 3) 5 ab + n a n b 



p a = E a (1 + — — W ) - — G E / £ b (T, X') v ' \J v ;, d X' + <D(c~ 4 ) , (33) 



2c 2 Jv v ' ' IX -X 



"" = CT (34) 

the PPN inertial and tidal torques read 
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inertial = J y X b (s Q c + ^ 2(1 + 7 ) £ C(fe (c d £ £ + 1 £ C d X e ) ) d 3 X, (35) 

^ = £flk |x 6 ^ TiC + i2(l + 7 ) ^ E A W c + sd (w^ d -W^ c )Jj d'X + Oic' 4 ), (36) 
and the additional torque L^ OT is defined as 

£n ot = e abc J y SX 6 f c rf 3 l = ^ e abc Q b E a% + 0(c~ 4 ) , (37) 

tt E = -- [ ZW E X a d 3 X + 0{c~ 2 ). (38) 
2 Jv 

The inertial torque is due to inertial forces dependent on the choice of Q a and C a . The 
PPN tidal torque is due to external gravitational potentials and It vanishes 
for an isolated body and is proportional to the square of the spatial size of the body. The 
additional torque L^ or represents an analogy to the Nordtvedt effect (known from the PPN 
translational equations of motion of extended self-gravitating bodies) in the rotational equa- 
tions of motion. L^ or is proportional to the Nordtvedt parameter rj, which is not zero in a 
particular theory of gravitation only if that theory leads to the violation of the Equivalence 
Principle. The effect results from the difference between the center of inertial mass and 
that of gravitational mass (the latter involves also the gravitational binding energy). The 
Nordtvedt effect consists actually in the fact that the two forms of mass-energy experience 
different accelerations provided that the Strong Equivalence Principle is violated. This fact 
results in the additional torque L^ OT . This additional torque in the rotational equations of 
motion enables us in principle to test the Strong Equivalence Principle with observations of 
rotational motion of celestial bodies: if Af3 — 7 — 3 7^ rotational motion of a body depends 
on its acceleration relative to the global reference system. L^ or is proportional also to the in- 
tegral n% which is an analogy to the internal gravitational energy of the body Qe appearing 
in the PPN translational equations of motion. Q E obviously vanishes for spherically sym- 
metric bodies. Roughly speaking, for Q E to be non-zero the body should be non-symmetric 
with respect to its center of mass (for example, 0| ^ for a homogeneous semi-sphere). 
This implies that £l E is quite small for typical celestial bodies. For that reason it is un- 
clear if this "new effect" appearing in the metric theories of gravity will lead to measurable 
consequences. 

More details about the local PPN reference system as well as about various equations of 
motion and their multipole expansions will be published elsewhere. 
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